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ABSTRACT 


'^'^Several  new  and  improved  shallow-water  normal-mode  models  have 
been  developed  to  represent  the  shallow-water  area  encountered  on 
the  FASOR  I and  II  criiises.  These  models  include  a ''’linear'* 
gradient  velocity  profile,  an  Epstein  velocity  profile,  and  a 
profile  with  three  Minear**  gradient  segments.  In  all  cases  a 
layered  viscoelastic  bottom  is  used.  Generally  good  agreement  is 
shown  between  the  calculated  and  experimental  m.easurements  of  propa- 
gation loss  on  the  FASOR  stations. This  memorandum,  has  been  prepared 
because  it  is  believed  that  the  information  may  be  useful  in  this 
form  to  others  at  KUV/C  and  to  a few  persons  outside  MJWC.  This 
memorandum  should  not  be  construed  as  a report  since  its  function 
is  to  present  information  that  later  will  be  3upplem.ented  and  ex- 
panded in  a report. 
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INTROLUCTION* 

For  shallow  water  propagation  at  noderate  sea  states  the  n-ost 
in.portar.t  environ  men  ta]  factors  ere  the  nature  of  the  bottom  and  the 
sound  velocity  profile.  At  KUVi'C,  San  Diego,  we  have  worked  to  develop 
models  that  accurately  represent  the  profile  and  the  bottom  sediments. 

We  will  use  these  models  for  calcvilauicns  of  propagation  in  certain 
areas  of  interest,  as  checks  for  approximate  shallow  water  models,  and 
to  study  special  aspects  of  propagation  such  as  the  effect  of  the  thermo- 
cline. 

We  will  present  two  models.  In  one  model  there  are  either  one  or 
three  water  layers  where  the  square  of  the  index  of  refraction  is  a 
linear  function  of  depth.  Since  the  velocity  is  almost  a linear 
function  with  depth  we  will  call  this  a linear  gradient  model.  This 
type  of  profile  has  been  used  by  Marsh^,  Tolstoy  , Pedersen  , and  others 
in  other  normal  mode  calculations.  The  second  model  is  a one  layer 
Epstein  profile.  Since  general  normal  mode  theory-  is  in  the  literature 
we  will  only  touch  on  the  unique  features  of  the  present  models. 


*This  technical  note  includes  part  of  the  information  presented  as 
Paper  G-8,  25th  U.  S.  Navy  Symposium  on  Underwater  Acoustics, 
Orlando,  Florida,  November  7-9,  1967.  This  work  was  accomplished 
under  NAVSHIPS  Sub-project  SF  101  03  15  Task  8105. 
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LINEAR  GRA.DIEtri'  I-.CLEL 


Figure  1 shows  the  potential  functions  for  a single  linear  gradient 
model. 

The  functions  h]_  and  h2  fcre  rodified  Henkel  fionctions  of  order 
one-third.^  The  bottom  is  made  up  of  an  arbitrary'  number  of  viscoelastic 
layers.  In  order  to  use  an  old  computer  subroutine  that  calculates  the 
reflection  coefficient  for  this  type  of  bottom.,  a pseudo  constant-velocity 
water  layer  is  set  in  between  the  depths  designated  z^'  and  z^,*  After  the 
usual  interface  conditions  are  set,  we  let  z^,' approach  zj-  so  that  the 
iso-velocitj'  layer  disappears  from,  the  problem.. 

The  argiur.ent  of  the  modified  Hankel  function,?;  > is  a function  of 
the  node  wave  number  k and  z.  The  perar.eters  a and  R are  functions  of  k. 

Finding  the  modes  is  the  mathematical  problem  determining  that  set 
of  k's  w'hich  satisfy  all  boiumdarj’  conditions.  For  example,  we  can  use 
the  boundary  conditions  at  z = z^  to  fix  A and  E.  The  mode  condition, 
also  called  the  dispersion  equation,  can  then  be  written  as  the  condition 
that  the  potential  he  zero  at  the  surface,  z = C.  The  dispersion 
eqmtion  is  stated  in  Figure  2. 

The  asymptotic  form  of  the  range  dependent  part  of  the  potential 
function  contains  ohe  term.  Thus  the  real  part  of  k determines 

the  horizontal  phase  velocity  of  a mode  and  the  imaginary  part  of  k 
fixes  the  attenuation.  n denotes  the  mode  nu-mber.  If  we  equate  the  real 
part  of.k  with  w/c',  where  c / is  the  horizontal  phase  velocity  of  a ray 
in  the  channel,  we  can  make  a connection  between  ray  and  mode  theory 
that  will  prove  helpful  in  the  solution  and  in  the  interpretation  of 


the  results.  The  value  of  Cj*  is  set  by  requiring  the  ray  to  satisfy 
a constructive  interference  relation ^ as  it  travels  down  the  channel. 

As  shown  in  Figure  2,  X is  the  horizontal  cycle  distance,  T is  the  tire 
j to  cor.plete  a cycle,  u is  the  angular  frequency,  is  the  phase  shift 

on  bottom  reflection,  €3  is  the  -^2  ® turn-over,  or  -/7  if  the 

ray  reflects  from  the  surface,  and  n is  the  rode  nurher.  The  approxi- 
mate value  of  the  iraginary  part  of  the  phase  velocity  is  calculated 
by  noting  the  loss  per  bottom  reflection  then  distributing  this  loss 
as  an  exponential  function  of  r. 

The  solutions  of  the  exact  and  approximate  dispersion  equation 
are  shown  in  Figure  3.  V»'e  have  plotted  values  of  node  velocity  c^rwAn 
rather  than  the  so  that  the  real  part  of  Cj-^  can  be  associated  with 
the  horizontal  phase  velocities  of  the  associated  rays  as  shown  in 
Figure  2.  For  example,  if  the  real  part  of  c^,  is  close  to  the  velocity 
at  the  bottom*  of  the  channel  the  node  will  be  confined  to  depths 
where  the  sound  velocity  is  less  than  the  real  part  of  Cj^.  When  the 
real  part  of  c^^  is  equal  to  the  surface  velocity  the  moce  begins  to 
reflect  from,  the  surface  and  the  attenuation  goes  up  because  the  cycle 
distance  X begins  to  diminish  in  length. 

The  open  circles  in  Figure  3 are  the  approximate  values  of  phase 
velocity  and  the  crosses  are  exact  values.  The  two  sets  of  Cjj  are 
essentially  identical  except  for  three  of  the  modes.  Note  that  the 
■ first  15  m.caes,  those  which  don't  reflect  from  the  surface,  have  about 
the  sam.e  attenuation.  This  is  because  as  the  .mode  number  goes  up  the 
bottom,  loss  goes  up  but  the  distance  between  reflections  also  becom.es 
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greater.  This  effect  of  loss  independent  of  grazing  angle  vas  used  by 

7 

Urick  in  his  shallow  water  ray  theory  rodel.  Later,  results  will  be 
discussed  for  a 3 layer  linear  gradient  case.  The  theory  is  essentially 
the  same  as  for  the  single  layer  model. 

EPSTEIN  MODEL 

The  second  model  concerns  a one  layer  Epstein  profile  as  shown  in 
Figure  4.  The  case  here  is  a transition  layer  with  the  same  surface 
and  bottom,  velocities  as  in  the  linear  gradient  mocel.  The  profile 
nay  also  contain  a hyperbolic  secant  term.,  however  this  case  will  not 
be  discussed. 

The  solution  is  similar  to  the  linear  gradient  model  except  the 

I potential  fvjictior.  is  new  a product  of  an  exponential  term,  tim.es  a 

i 

; hypergeonetric  function.  The  exact  form,  of  these  term.s  can  be  fc  nd  in 

e 

a paper  by  Deavenport  . 

I 

Values  of  phase  velocity  for  this  profile  are  shown  in  Figure  5. 

In  this  case  the  mode  attenuation  increases  as  the  boutor  loss  in- 
creases but  then  has  a sharp  decrease  as  the  phase  velocity  approaches 
the  sound  velocity  in  the  upper  part  of  the  channel.  The  result  is  two 
sets  of  low  attenuation  modes.  The  first  group  strong  below  the  thermo- 
cline  and  the  second  set  strong  above.  Thus  propagation  across  the 
thermocline  suffers.  For  the  present  values  of  surface  and  bottomi 
velocity  and  the  case  of  an  extremely  sharp  ther::ocline  (two  iso-velocity 
layers)  the  loss  across  the  thermocline  is  5 to  10  db.  This  is  at 
i l.JkH^,  t’.  .■  loss  woiild  be  greater  at  higher  frequencies.  The  rule 

' is:  don't  cross  your  moc.es  in  the  rlcdle  of  a propagation  run. 


MODE  3U>i-ATI0N 


To  calcultte  the  scund  field  the  nodes  ere  sunr.ed  as  shoun  In 
Figure  6.  The  second  forin  of  the  norir.alizing  term  can  generally  be 
calculated  vd.th  better  accuracy  because  it  is  depender.t  upon  ^ calcu- 
lated throughout  the  channel,  rather  than  only  at  the  surface.  To 
facilitate  the  integration  of  over  the  sediinent  layers  we  use  a trick 
as  shown  in  Figure  7. 

Taking  the  nodes  one  at  a time,  we  pick  an  absorbing  liquid 
half-space  that  has  the  same  density  as  the  water  and  will  give  the 
same  reflection  coefficient.  That  is,  the  layered-viscoelastic  bottom 
is  used  to  find  the  modes  but  we  switch  to  an  absorbing-liquid  half- 
space for  the  normalization.  The  integral  from  z^  to  oo  can  be 
completed  by  inspection  since  there  is  only  an  exponential  wave  in  the 
absorbing  liquid. 

RhoULTi  - 3 LtYlR  LIM-aR  GRr.LIEl.'T  KOLhL 

There  were  three  shallow  water  stations  on  the  last  Fa3CR  cruise 
in  the  Pacific  by  I.UV.'C,  Cn  one  of  these  the  tv.o  isc-veiocity  layer 
model  provided  a good  fit  to  the  velocity  profile.  On  the  other  two 
stations,  OAK  and  THCRIs,  we  have  made  intensity  calculations  using  a 
three  layer  linear  gradient  model.  In  Figure  8 are  ET's  taken  at  the 
tim.e  of  the  tests  for  station  OaK.  The  tests  extended  over  17  hours. 
The  broader,  heavier  line  is  our  fit  to  the  profile. 

Measured  and  calculated  vr lues  of  propaga mien  loss  for  ranges 
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up  to  50  kyd  are  shown  in  Figure  9.  Thorpe's  value  of  volum.e  attenua- 
tion and  Amos  values  of  surface  loss  for  sea  state  1 were  added. 
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That  is,  we  calculated  propagation  loss  from  Karsh  and  Shulkin's 
equations  for  sea  state  1 and  sea  state  0.  The  difference  was  taken 
as  the  surface  loss.  The  agreement  between  calculated  and  experimental 
results  is  somewhat  better  than  we  reported  last  year  using  a two  layer 
iso-velocity  rcdel.  However,  the  experimental  losses  are  about  6 db 
greater  over  the  section  of  the  range  interval  from  2C  to  4C  kyd. 

The  BT's  and  the  three  layer  .lOGel  for  THORN  are  shown  in  Figure  10. 
The  profile  is  interesting  because  of  the  well  developed  surface  duct. 
The  modes  separate  namurally  into  three  groups.  A set  that  corres- 
ponds to  rays  that  have  only  bottom  reflections,  two  roces  in  the 
surface  duct,  and  a group  of  higher  order  modes  that  reflect  alternately 
from  the  surface  and  bottom.  The  mode  phase  velocities  are  shown  in 
Figure  11.  Again  the  approximate  roots  are  an  accurate  estinate  of  the 
real  parts  of  the  set  of  c^-'s  for  all  modes  and  of  the  imaginary  parts 
for  all  nodes  but  one.  The  calculated  and  observed  values  of  propaga- 
tion loss,  Figure  12,  can  be  seen  to  be  in  good  agreement. 

CONCLL'SIOM 

k'ith  the  present  models  we  feel  that  a satisfactoim’  first  order 
theory-  has  been  developed  for  shallow  water  propagation,  Cn  this 
fram.ework  additions  can  be  added  to  account  for  boundary  roughness, 
horizontal  variations,  and  other  second  order  effects. 
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EXACT  DISPERSION  EQ-: 

4>,(z  = 0,  k = kn)  = 0,  n=  l,2,3,--- 

APPROXIMATE  DISPERSION  EQ.: 
kn=  w/cn“w/(cn-icn  ) 

w(T“  X/Cp)  + + €g  = (n-l)2ir 

Cfi  =.(Cf| ) • Bl / (x- 0) • 8-686) 


FIGURE 
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